Response and Amplification of Terahertz Electromagnetic Waves in Intrinsic Josephson Junctions 

of Layered High-r, Superconductor 



O 

(N 



Shi-Zeng Lin and Xiao Hu 
' WPI Center for Materials Nanoarchitectonics, National Institute for Materials Science, Tsukuba 305-0044, Japan 
-Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan 

(Dated: July 13, 2010) 

We investigate the response of a stack of intrinsic Josephson junctions (IJJs) to terahertz (THz) electromag- 
netic (EM) irradiation. A significant amplification of the EM wave can be achieved by the IJJs stack when the 
incident frequency equals to one of the cavity frequencies. The irradiation excites n phase kinks in the junc- 
tions, which stimulate the cavity resonance when the bias voltage is tuned. A large amount of dc energy is then 
pumped into the Josephson plasma oscillation, and the incident wave gets amplified. From the profound current 
step in IV characteristics induced at the cavity resonance, the system can also be used for detection of the THz 
wave. 

PACS numbers: 74.50.+r, 74.25.Gz, 85.25.Cp 
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It has been known for a long time that Josephson junc- 
tions can be used as oscillator, amplifier and detector for 
electromagnetic (EM) waveHl]. The operating frequency of 
these devices made of conventional low-temperature super- 
conductors is below terahertz (THz) due to the small super- 
conducting energy gap. The discovery of intrinsic Joseph- 
son effect in layered high-Tc superconductors [2||, such as 
Bi2Sr2CaCu208+5(BSCCO), has extended the frequency to 
the THz band, where EM waves have potential for wide 
applicationsisll^l, and thus has stimulated intensive research 
activities in the field fsl]. 

A breakthrough in generating coherent THz emission has 
been achieved recently based on a mesa structure of BSCCO 
single crystal [6]. Due to the thickness much smaller than the 
wave length of EM wave, the IJJs stack itself forms a cavity, 
synchronizes the plasma oscillation and radiates coherent THz 
wave at the cavity resonance ff\. The dynamics of the super- 
conductivity phase has been addressed theoretically that 
+n phase kinks are developed in the junctions, which couple 
the dc bias to the standing wave and make the cavity resonance 
possible. 

In applications THz detector and amplifier are as important 
as generator. Shapiro steps ifioll were observed in IJJs stacks 



of small BSCCO mesas under THz irradiation 111 1141411 . which 
can be used for detection. To the best of our knowledge, no 
experiment on amplification of THz waves based on IJJs has 
been reported so far. Now with the success of generation of 
coherent THz wave at cavity resonance it is intriguing to 
explore the possibility of the same setup for the usage of am- 
plification and detection of THz waves, with the expectation 
that the system exhibiting a cavity resonance in the THz band 
responds more sensitively to an incident wave than short junc- 
tions reported in literatures. 

When a stack of IJJs is irradiated by an EM wave, the trans- 
mitted wave excites Josephson plasma oscillations inside the 
IJJs. The incident wave can be either damped or amplified 
according to the detailed compensation between dissipations 
caused by quasiparticles and the power supply from the bias 
voltage, which, in turn, is governed by the phase dynamics in 



the stack of IJJs. 

By investigating the inductively coupled sine-Gordon equa- 
tions under appropriate boundary condition taking into ac- 
count the THz EM irradiation, we show in the present article 
that with the IJJs stack one can achieve a significant amplifi- 
cation of the input wave with frequency equal to the one of the 
cavity frequencies. Tuning the bias voltage, n phase kinks are 
created in the junctions, which pumps a large amount of dc 
energy into the Josephson plasma oscillation due to the cav- 
ity resonance. The profound cuiTent step in IV characteristics 
induced at the cavity resonance signals the existence incident 
THz wave, and thus can be used for detection. 

The setup is shown in Fig.[Tl where a stack of IJJs are sand- 
wiched by two ideal conductors with infinite thickness. These 
two conductors prevent the interference between EM waves 
from the two edges of IJJs stack. The left side of IJJs is ex- 
posed to irradiation. We assume that the IJJs are infinitely 
long in the y direction, and thus the problem reduces to two di- 
mensions with sizes and Lj. This setup is similar to the one 
proposed in Ref.fisl, except for the lateral size - lOOyum 
which contains cavity modes in the THz regime. 

The dynamics of the gauge invariant phase difference in 
IJJs is described by the inductively coupled sine-Gordon 
equations||l[i[il[lJ 



dlP, = (1 - ^Ad)[sinP, +l3d,P, + djP, - /ex,], 



(1) 



where Pi is the gauge invariant phase difference at the l-th 
junction, /3 = Ano-cAe/c y/s^ the normalized c-axis conductiv- 
ity, ^ - (Aah/s)^ the inductive coupling; is the dielectric 
constant and cr^, is the conductivity along the c-axis, and s 
is the lattice period in the c-direction; c is the light velocity 
in vacuum; Ac and Agt, are the penetration depths along the 
ab-axis and c-axis respectively. In Eq.O, the lateral space 
is normalized by A^, time by the Josephson plasma frequency 
(jjy = c/Ac yfe^, and the external current 7ext by the Josephson 
critical current 7c llSll . Ad is the second-order difference op- 
erator defined as Aa/z = fi+\ -n fi-\ - 2//. We adopt fi - 0.02 
and ^ = 7.1 X 10^ which are typical for BSCCO [19]. The 
physics discussed below is valid in a stack of IJJs with huge ^, 
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FIG. 1; (Color online) A stack of IJJs sandwiched by two electrodes 
and subjected to irradiation. The whole setup is immersed in a di- 
electric material with dielectric constant sj. 



which has not yet been achieved in artificial Josephson junc- 
tion stacksl 16i 20ll. 

In the absence of irradiation, the so-called dynamic bound- 
ary condition (DBC) was derived based on Maxwell equations 
10 lilt]. It is easy to generalize the DBC to incorporate irra- 
diation because Maxwell equations are linear. Assuming the 
incident wave is a plane wave with the electric field polarized 
along the z axis and the propagation direction normal to the 
left edge of IJJs, the total electric field on the left side is 



E^ix, z, t) - E\{lli) exp[/(- -\/eda)x + tot + 8)] 
+ J dkxE°{u, kx, k,) e\p[i(-k^x + k,z + cjt)]. 



(2) 



where k\ + - oj^ with E° the outgoing wave comprising 
the emitted and reflected waves, E\ the incident wave with a 
relative phase difference Q to the Josephson plasma oscillation 
inside IJJs, w the frequency, and the normalized dielectric 
constant of the dielectric medium coupled to the IJJs ifisll . For 
simplicity of analysis, we concentrate in Eq.© on the case 
that the frequency of Josephson plasma determined by the bias 
voltage according to the ac Josephson relation is equal to the 
incident frequency, since otherwise the response of IJJs is very 
small. 

The EM wave at the right edge comes only from emission. 
The generalized boundary conditions for the oscillating elec- 
tromagnetic fields in the real space and frequency domain are 
given bylfisl] 



By(X = 0, uS) 



EAO, oj) 
Z(w) 



2 V^£j(0, a>) expiiO), (3) 



By{X = L^, u) 



E-SL„oj) 
Z(w) ' 



(4) 



where By and E- are the total electric and magnetic fields, 
and Z(w) = 2/ {fe„L, V^[l + f ln||]} with k^ = ujyf^ 



is the impedance 10] • The power of the incident wave is 
5'i = y[ed(_E\)^ 12. Since the thickness of the IJJs stack used in 
experiments is L, - lyum, much smaller than the wave length 
in the THz band, the electromagnetic fields are uniform in the 
z direction in the IJJs for THz waves as in Eqs. Q and (|4|. 
There exists a significant impedance mismatch between the 
IJJs and the dielectric medium, which is crucial for the cav- 
ity resonances. With the relations (1 - fAd)B) - d^Pi and 
one obtains the boundary condition for the os- 
cillating part of Pi. 

A solution to Eq.O is given intuitively by 



Pi{x, t) - (jjt + Re[-ig(x) exp(/ajO], 



(5) 



where the first term at the right-hand side (r.h.s) is the uniform 
rotating phase according to the ac Josephson relation and the 
second term is the plasma oscillation. The spatial modulation 
of plasma oscillation is induced by both radiation and irra- 
diation. We consider the region of small plasma oscillation 
\g{x)\ < 1. From terms with time dependence exp(+i(L)t), we 
obtain the equation foig{x) by substituting Eq.Q into Eq. ([T]l 



dlgix) = 1 + i/3ajg(x) - g{x)a?. 
Equation (|6]l has the solution 

g(x) - A + a exp(iqx) + b exp{-iqx) 



(6) 



(7) 



with A = 1 / {or - ijiuS) and q ■x. cu for weak damping yS <«; 1 
as in the case of BSCCO system. The first term A in Eq. (|7J 
represents the uniform plasma oscillation, and the other two 
terms are propagating waves due to radiation and irradiation, 
with the two coeflicients a and b determined by the boundary 
condition Eqs. (|3]l and (|4|i. 

The TV characteristics is derived from the current conserva- 
tion relation 



:^w + <sinP,>^, 



(8) 



where (■ ■ ■ )^., denotes the average over space and time. Be- 
sides the normal current due to the quasiparticles - 
the total dc current has two contributions from the plasma 
oscillation due to the nonlinearity of the dc Josephson ef- 
fect: 7p = I3l[2(a? + 13^0))] and J„ = {a[\ - exp(/wLi.)] - 
b{\ - exp{-i(jL>Lx)\] I (2(jL>Lx) associated with the uniform and 
nonuniform parts of plasma oscillation in Eq.©. The IV 
curves for different 0's are displayed in Fig. |2a) for the in- 
cident wave 5'i = 141W/cm^. When one fixes the voltage 
satisfying the phase-locking relation and sweeps the current, 
the relative phase adjusts itself to match the current, which 
traces out the Shapiro steps [IQ,]. Zero-crossing Shapiro steps 
II22II occur at small voltages. For 1/Z « coL,^ 1, is given 
explicitly as 



- Re 



1 EW^ 



(9) 



/w is maximized (minimized) ?A6 - ji {9 - 0), and the height 
of the Shapiro step is given by - 2E\ yfedlL^a?. In the 
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FIG. 2: (Color online), (a) W characteristics and (b) radiation power 
at tiie right edge of the state described by Eq. ([Sj for the incident 
wave of 5, = 141W/cm^ at several typical phases Q. The light yellow 
regime in between the maximal and minimal current associated with 
= and Q = n '\^ the width of the first Shapiro step. The results are 
obtained with = 80/jm, L- = Ifim and = 0.1 similar to those in 
the experiments 



region where AJ^ = Jv^iEl > 0) - /„(£'' = 0) > 0, a dc power 
A7w<j^ is converted into emission. When 7w < 0, the incident 
EM wave is converted into dc power and charges the IJJs, and 
the IJJs effectively work as a battery. 

The radiation powers measured by the Poynting vector sat- 
isfy the power balance condition S r - S i + - /ext'^ with P^i 
the dissipation caused by the quasiparticles, and Si andS r the 
radiation at the left and right edge respectively. The emission 
at the right edge is depicted in Fig.|2b) for Si - 141W/cm^. 
It is clear that the emission is always very weak because of 
lack of an efficient way to pump energy into plasma oscilla- 
tion in this state. In the same limit, the radiation power is 
given by 



Sr = Re 



2Z* 



2iE\e"' ^|^d 



(10) 



It is clear that the emission comprises of spontaneous one 
•Ssp ~ I^P, the one caused by transmitted wave S'tw ~ \E['^ 
and the stimulated one ^st ~ \E'^\. 

The state given in Eq.© is found to become unstable when 
the incident frequency and the bias voltage get close to u - 
k\ = n/Lx, where the cavity mode g(x) x A[ cos(^i.x) with 
Ai > 1 is induced by the irradiation (without losing generality, 
here we consider the first cavity mode). Instead of Eq.© the 
phase dynamics should then be described by lit] 

Pi(x, t) - Lot + P'i(x) + Re[-ig(x) exp(/<yf)], (1 1) 

where P^j is the static phase kink associated with the cavity 
mode. Substituting Eq. (fTTT i into Eq. ([T]i, we obtain the fol- 
lowing equations 

dlg(x) = (1 - f Ad) exp(iPj) + ipcjgix) - g{x)u? (12) 



FIG. 3: (Color online), (a) W characteristics and (b) radiation power 
when the n phase kinks are excited near the cavity resonance. Other 
parameters are the same as those in Fig. |2] 



and 



= -fAdg(^)exp(-/P,0. 



Since the Josephson plasma should take the form 
g{x) - A\ cos k\x + a QYjp(iqx) + h exp(-/^x), 



(13) 



(14) 



where q x, cu, ys/e. arrive at Ai = Fi/{ikf - ia> - /3(d) with 
Fi - 7^/(1 - fAd)exp(/Ppcos(fciJc)iix. Due to that the 



dominant term in Eq.(fT4l) g{x) x Ai cos(feix) is antisymmet- 
ric with respect to x - L J2, Eq. ( fT3l l has a solution of n 
kinks alternatingly piled up in the c-axis|8]. Because of the 
huge inductive coupling ^ ^ 10^ in BSCCO, the phase kinks 
render themselves as step functions, and are stable against the 
radiation and irradiation provided |fl|, |/7| < |Ai|. The width ofa 
kink should be smaller than the junction width 1 / -y/|Ai|^ < L i , 
which gives an estimate on the regime where the kink state is 
stable. 

Although it has been revealed theoretically that the n kink 
state is ideal for generating strong terahertz electromagnetic 
waves [5, 8], the dynamic process to realize the state was not 
clear. The present study indicates that irradiating the junction 
stack by an incident wave can stimulate the tt kink state. 

The dc supercurrent induced by the plasma oscillation in 
the kink state can be evaluated by (^-ig(x)e\p(-iPJ)/2^ . It 
includes the current associated with the plasma oscillations at 
the cavity mode Jp - A-jSo) I [n\{k\ - ap-)^ + 0^(jJ^]], and that 
with radiation and irradiation - (a [exp(/dLtLv/2) - 1]^ - 
b [exp(-/wLt/2) - lf')l{2LoLx). The IV characteristics with 
the irradiation of Si - 141W/cm^ is given in Fig.^a). For 
1 /Z <s; A:i - w «: 1 and {k\ - a?) » pco, we have 



7w = Re 



2E}/^ 



L.n^ZiB - w2)(fei - (jj) n\ky - oS) 



(15) 
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FIG. 4: (Color online). Amplification factor for several typical inci- 
dent powers near the cavity resonance. 



Because of the n phase kinks, now 7„ is maximized (mini- 
mized) 6 - Q (6 - 7i) m contrast to Eq. (|9j. The height 
of the Shapiro step is 7s = 4£'' yfe^lin^ikx - w)]. As is well 
known, Shapiro steps are suppressed by internal modes in sin- 
gle junctions. It is the same case for a IJJs stack if the state is 
uniform along the c-axis, since Eqs.([T]i are decoupled. There- 
fore, the appearance of a Shapiro step at the cavity resonance 
can be used as an exclusive detection of the n kink state in a 
stack of IJJs. 

The radiation power at the right edge is depicted in Fig. [3b) 
for Si - 141W/cm^. It is clear that the input wave is en- 
hanced significantly near the cavity frequency, where the n 
kinks stimulated by the irradiation pump a large amount of dc 
power from the dc bias into Josephson plasma oscillation. 

The radiation power at the right edge in the same limit is 
given by 



Sr = Re 



2Z* 



7l{k] - W^) T^ihx - 



(16) 



An amplification factor can be defined by the maximal value 
of the ratio S r/Si with respect to the phase 6 for a given Si. 
For the tt kink state, one has 



^/edZ* 



2 ^/sd 



(17) 



which is achieved al9 - njl. As displayed in Fig.|4] the am- 
plification factor reaches its maximum at w = An incident 
wave of power of 141 W/cm- can be amplified by one order of 
magnitude. The amplification factor decreases with the power 
of incident wave, and the maximum power which can be am- 
plified by this technique is estimated as 3000W/cm^. 

In a single junction, the presence of irradiation may cause 
chaotic dynamics in a certain parameter space ll22ll . which is 
harmful for applications. The chaos can be avoided when the 
frequency of the incident wave is much larger than the Joseph- 
son plasma frequency] 12 22 1, which is fulfilled in a IJJs stack 
of length smaller than Ac- 



In conclusion, simultaneously shining a terahertz electro- 
magnetic wave and biasing a dc voltage on a stack of intrinsic 
Josephson junctions stimulates the standing wave of Joseph- 
son plasma, which develops n phase kinks in the junctions, 
when the frequency equals to one of the cavity frequencies 
of the junction stack. At the cavity resonance, the rotating 
n kinks pump a large amount of dc energy into Josephson 
plasma oscillation, and the incident wave gets amplified. The 
maximal radiation power reached by this terahertz amplifier 
is estimated as 3000W/cm". Since the strong plasma oscilla- 
tion induces a large dc supercurrent at the cavity resonance, 
the system can work as a terahertz detector The response of 
the system to irradiation depends on the spatial structure of 
the superconductivity phase, thus the phase dynamics may be 
probed by the irradiation. 
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